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ABSTRACT 

The meaningful comparison of models of galaxy evolution to observations is critically dependent on the ac- 
curate treatment of dust attenuation. To investigate dust absorption and emission in galaxies we have assembled 
a sample of ^1000 galaxies with ultraviolet (UV) through infrared (IR) photometry from GALEX, SDSS, and 
Spitzer and optical spectroscopy from SDSS. The ratio of IR to UV emission (IRX) is used to constrain the 
dust attenuation in galaxies. We use the 4000A break as a robust and useful, although coarse, indicator of star 
formation history (SFH). We examine the relationship between IRX and the UV spectral slope (a common at- 
tenuation indicator at high-redshift) and find little dependence of the scatter on D„(4000). We construct average 
UV through far-IR spectral energy distributions (SEDs) for different ranges of IRX, D„(4000), and stellar mass 
(M„) to show the variation of the entire SED with these parameters. When binned simultaneously by IRX, 
D„(4000), and M* these SEDs allow us to determine a low resolution average attenuation curve for different 
ranges of ■ The attenuation curves thus derived are consistent with a A attenuation law, and we find no 
significant variations with M* . Finally, we show the relationship between IRX and the global stellar mass sur- 
face density and gas-phase-metallicity. Among star forming galaxies we find a strong correlation between IRX 
and stellar mass surface density, even at constant metallicity, a result that is closely linked to the well-known 
correlation between IRX and star-formation rate. 

Subject headings: galaxies:fundamental parameters — galaxies:evolution — dust:extinction — ultravio- 
let:galaxies — infrared:galaxies 



1. INTRODUCTION 

Measurement of the star formation history (SFH) of galax- 
ies, and the distribution thereof, tests models of galaxy evo- 
lution. Such measurement requires the detailed treatment of 
dust attenuation, which poses one of the most severe obsta- 
cles to converting observed properties (e.g. broadband colors 
or line fluxes) into a SFH or star formation rate (SFR). 

One of the primary tools for the measurement of dust at- 
tenuation, especially at high-redshift, is the relation between 
the ratio of infrared flux to ultraviolet flux (IRX) and the ul- 
traviolet spectral slope (fi, where f\ ~ A' 3 ). The former can 
be considered an approximate measure of the amount of dust 
attenuation since the IR flux is produced by light absorbed at 
primarily UV wavelengths, while the UV flux measures the 
light transmitted at these wavelengths. The UV spectral slope 
measures the reddening (or color excess) of the UV spectrum 
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due to selective absorption by dust - assuming that before at- 
tenuation the spectrum is nearly flat and relatively insensitive 
to the SFH. The utility of this relation lies in the relatively 
easy measurement of the restframe UV color for large sam- 
ples of galaxies at z > 1, from which the attenuation can be 
inferred. 

The calibration (and scatter) of this relation is of great inter- 
est for measurements of the SFR of galaxies at high redshift 
(especially those based on the restframe UV luminosity) and 
for t he resulting estima t es of the SFR density and its evolu- 
tion dMadau et al.ll 19961: [Schiminovich et alj|2005l) . This cal- 
ibration is usually empirical, as is necessary at least in part 
because of the dependence of the reddening effect of the dust 
on the assumed properties of the dust (e.g., the extinction law 
and dust geometry). One must then assume that these dust 
properties do not evolve strongly with redshift, although cal- 
ibrations based e ntirely on high-red shift galaxies may soon 
become available (Reddy et al. 2006|). 

Such empirical calibration, at any redshift, requires both 
restframe IR and UV observations galaxies. In the past, such 
data has been available from IRAS and IUE for samples of 
nearby starbursting galaxies, w here a strong correla tion be- 
tween IRX and (3 was found dMeurer et all 1 19991) . With 
the launch of GALEX, UV data has bec ome available for a 
larger sample of more normal galaxies (ISeibert et all 120051 : 
iBuat et all 120051) . for which it appears that the relation be- 
tween IRX and (3 is shifted to redder colors, and has increased 
scatter. Similar results were obtain ed for a sample of individ- 
ual OB associations by lBelll (H002), who suggested that these 
effects we re due to variations in the relative star-dust geom- 
etry (e.g., IWitt & Gordonll2000h . iKong et all (12004 found a 
trend with the 4000A break (a coarse measure of the SFH) of 
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the offset of starburst galaxies from the best fit IRX-/3 relation. 
On the basis of spectral synthesis modelling they thus ascribed 
the scatter and shifts in towards redder colors in the IRX-/3 re- 
lation to the effect of SFH on (3, in the sense that a large popu- 
lation of old stars will tend to produce a redder intrinsic, unat- 
tenuated, UV color. The largest effects in the models were 
seen for significant old bursts superimposed on a smoother 
SFH. The spatial resolution of GALEX, combined with the 
spatial resolution of Spitzer, has allowed for the investiga- 
tion ofjhe_TRX^rehyionj^ 

ies dBoissier et alj|2004: iGil de Paz et alj|2007t ICalzetti et all 
120051) . One important outcome of this work is the sugges- 
tion that populations of different ages ar e obscured by differ- 
ing amounts of dust (Calzetti et al. 2005), a result consistent 
with the findings o flCharlot & F all (2000) for global measures 
of attenuation. Such a scenario f urther complicates the inter- 
pretation of the IRX-/3 relation dPanuzzo et all 120071) . The 
combination of sensitive UV and IR observations of galaxies 
with GALEX and Spitzer has enabled the in vestigation of the 
IRX-/3 relation for l arge samples of galaxies (Buat et al .120051 : 
ICortese et al]|2006l) . 

This paper is part of a series investigating the UV through 
IR properties of a large, well-defined sample of galaxies 
observed spectroscopically for SDSS, and by GALEX and 
Spitzer. The additional diagnostics of stellar populations, 
SFR, and attenuation provided by the SDSS data make the 
sample presented here a unique and important testbed for the 
understanding of the IRX-/3 relation. 

The large amount of homogenous ancillary information 
available from the SDSS observations, in addition to the UV 
and IR data, also allows us to investigate the effects of dust 
across the UV through near-IR spectrum. When such a wide 
range of wavelengths is considered, the effects of SFH on the 
spectrum become much more pronounced th an they are for 
the UV spectrum alone (John son" et al.ll2007a|). Indeed, opti- 
cal colors are ofte n used as a proxy for SFH (lBelleLaD[2004; 
iFaber et al.l2005l) . despite the potentially large additional con- 
tribution of attenuation to these colors. Here we present the 
global UV through IR spectral energy distributions (SEDs) of 
a large range of galaxy types. We show how these measured 
SEDs vary as a function of stellar mass, SFH, and attenua- 
tion. In particular, we use these SEDs to derive an 'average' 
dust attenuation law, and make a first attempt to determine 
the variation of this law (if any) with the stellar mass of the 
galaxy. These SEDs are complementary to the integrated UV 
through IR SEDs presented for the smaller s amples of much 
more local , resolved galaxie s presented by IGil de Paz et al.1 
(2007) and lDale et aD (12006b . and to the average UV th rough 
IR SEDs of higher redshift galaxies (IZheng et alj|2007l) . 

Besides constraining the correction of derived physical val- 
ues for dust attenuation, the relation of the dust attenuation 
to the SFH, SFR, and stellar mass may provide an additional 
constraint on models of galaxy evolution, as the attenuation 
should be proportional t o the gas surface density S g?t and 
gas-phase metallic i tv Z dWang & Heckmanl fl996: Belll l2003h 
iMartin et ail l2006t ICortese et al l 120061) . Models that self- 
consistently treat absorption and emission by dust are neces- 
sary and must be able to both predict the co rrect distribution of 
E(, aj Z (as are given by the models of, e.g.lCroton et alj|2006t 
iDe Lucia etaill2006t ISomerville et al.ll2001l) and convert this 
quantity into a dust attenuation curve to accurately predict the 
UV through IR properties of galaxies. We investigate the re- 
lation between T, gas Z and attenuation for the galaxies in our 
sample that have metallicity measurements. 



TABLE 1 
Observations 



Field Name 


Size 


N o4s 










(deg 2 ) 




f,n< 25 


m 2 4 < 19.5 




Lockman Hole 




872 


792 


819 


721 


FLS 


~3 


186 


147 


158 


118 



2. DATA 

In this study we use optical spectroscopic and photomet- 
ric observations of galaxies from SDSS, UV observations of 
these galaxies from GALEX, and IR observa tions of these 
galaxie s from Spitzer. Readers are referred to Johns onet alj 
(2007a) for details of the data reduction and sample defini- 
tion. Here we provide a summary of the sample properties 
and briefly describe several parameters derived from the data 
that are used throughout this analysis. 

The sample consists of galaxies targeted spectroscopically 
by SDSS (thus having a magnitud e r < 17.7), for which a 
stellar mass has been determined by Kauffm ann et alj d2003al) 
and which have been observed by both GALEX and Spitzer. 
This sample is located in two regions of the sky, the Lock- 
man Hole (observed by Spitzer as part of the SWIRE survey) 
and the Spitzer Extragalactic First Look Survey (FLS). Ta- 
ble Q] gives the number of galaxies observed in each field and 
the detection rates. The sample properties (redshift, stellar 
mass, and morp hological distribution, etc .) are very similar 
to the sample of Kauffm ann et al.1 d2003al) . although the red- 
shift distribution is less smooth due to the presence of groups, 
filaments, and sheets in the spatial distribution of galaxies in 
these fields. The median redshift of the sample is z = 0. 1 1 with 
a maximum of z < 0.305. The median stellar mass (see below) 
is logM* = 10.7M o with a range of 8. 45 < logM* < 11.91 . 
The SFR range is simO Ol-eOM^-yr- 1 dJohnson et alj|2007ri 
iBrinchmann et alj|2004l) . A detailed investigation of the mor- 
phological properties of the ample is beyond the scope of this 
work, but we note that irregular, elliptical, and spiral galaxies 
of v arious inclinations (fac e-on through edge-on) are included 
(see iJohnson et a l. 2007a, Fig. 4). The SDSS spect r oscopic 
galaxies with stellar masses from Kauffmann et al. ( 2003a) 
that are not detected in the UV or IR are ~ 30% of galax- 
ies on the red sequence and a few blue, low stellar mass dwarf 
galaxies that are missed primarily because we require objects 
to have Petrosian radii in the r band of less than 1 l"and red- 
shift z > 0.02, insuring accurate photometry. Photometry is 
from the modified SExtractor output for the GALEX far-UV 
(/, A - 1528A) and near-UV (n, A ~ 2271 A) data, SDSS 
Petrosian magnitudes (ugriz) for the optical data, and large 
aperture photometry for the Spitzer data (with 7"radius in 
the IRAC 3.6pm through 7.8/^m bands, 12"radius at 24pm, 
and 16"radius at 70pm where the sources are largely unre- 
solved). The / through 3.6/immagnitudes are /f-corrected to 
z = 0.1 (the median redsh ift of the sample) via the method of 
iBlanton & R oweis (2007), with the result denoted by a super- 
scripted 0.1, e.g. 01 m. We also determine Zf-corrections of the 
/ and n magnitudes t o z = 0.0 ( /, n) for comparison to 
more local galaxies ( j33.21 l, 

In this work we are in t erested in the stellar masses de- 
rived by Kauffmann et al. (2003a) and the metallicities de- 
rived by iTremonti et alJ d2004l) . These properties have been 
made available as catalogs 11 by these authors. The stellar 

1 1 http://www.mpa-garching. mpg.de/SDSS/ 
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masses are derived from the z-band SDSS photometry. The 
mass-to-light ratio in that band is obtained for each galaxy 
from fits of stellar population synthesis models to the opti- 
cal spectrum, in particular the 4000Abreak strength and the 
HS absorption strength, HSj^. The gas phase metallicities are 
derived from the SDSS measured emission lines for galaxies 
with strong emission lines that do not show evidence of AGN 

I III 

activity in the optical (see ITremonti et al.H2004L for details). 
We als o consider in §|6]the SFR derived by Brinchman n et al.l 
(2004) from a comparison of measured optical emission lines 
to a large suite of models of galaxy spectra. 

2.1. D„(4000) 

The 4000A break strength, here defined as in Balo gh et al.1 
(1999) (D„(4000)), is a useful indicator of the mass-to-light 
ratio since it is sensitive to the star-formation history (SFH). 
Larger D„(4000) signifies, approximately, a larger ratio of old 
stars to young stars. This break is measured at high signal- 
to-noise ratio (S/N) for all of the galaxies in the sample, and 
is not very sensitive to the presence of dust attenuation (but 
see lMacArthurf l2005 ) . D„(4000) is only measured within the 
3" radius spectroscopic aperture of the SDSS, and can there- 
fore b e overestimated for gala xies with moderate bulge/disk 
ratios (Kauffm ann et al.l 12007). D„(4000) serves as our pri- 
mary measure of SFH. 

2.2. IRX 

The combination of UV and IR data allows us to con- 
struct a robust measure of attenuation: the so-called infrared 
excess, IRX= log^^/L,,,,). We define L uv = vL v where 
v = c/1390A (A = 1390Ais the effective wavelength of the 
01 f band) and calculate L v from the /^-corrected absolute 
magnitude, assuming the redshift-distance relation given by 
a concordance cosmology with Q m = 0.3,0a = 0.7, Ho = 70 
km s" 1 Mpc" 1 . Ldus,, the 8-100 0/xm luminosity, is est imated 
from the 24/im luminosity as in Johns onet al.l (12007a). with 
a bolometric correction that depends on the ratio of 8/rni to 
24/im fluxes. These bolometric corrections are derived from 
the models of iDale et all (1200 ll) . 

If the IR emission traces the (predominantly blue) light 
from (predominantly) young stars that is absorbed by dust and 
reradiated, and the UV emission traces the light from young 
blue stars that is transmitted, then their ratio is a measure of 
the optical depth. The UV attenuation in magnitudes can be 
written in terms of IRX, including a parameter 77 to account 
for light absorbed by dust at wavelengths other than the UV, 
as 

A 1RX = 2.5log (r,10 IR X +l) (1) 

with r]= 1/1.68 dMeurer et alj|1999t Uohnson et al.ll2007bl) . 
Other conversions between IRX a nd attenuation are available 
(e.g.. lBellll2002l:lBuat et al.ll2005l) but are consistent with this 
formulation, which is physically motivated and simply under- 
stood. However, this formulation relies on a number of sim- 
plifying assumptions. First, the fraction of light absorbed by 
dust at wavelengths other than A ~ 1390A may depart from 
the canonical value of 1-1/1. 68, depending especi ally on the 
SFH of the galaxy (see, e.g.. lJohnson et al. 2007b) - this will 
change the relation between A/«x and the true attenuation at 
1390A, Ai39o- Related to this, different stellar populations 
withi n a galaxy may be attenuated by differe nt amounts of 
dust (ICharlot & Fallfl2000t ICalzetti et al1l2005l) . Second, the 
relative geometry of the stars and dust can lead to the mis- 
estimation of A 1390 by Amx, with an extreme example given 



by a group of s tars surrounded by a ring of dust normal to 
the line of sight ( Bel l et al.ll2002l) . This effect has been exam- 
ined in detail by IWitt & Gordon! d2000h. iBuat & Xul (fl996l). 
iGordon et all (120001) . and iPierinietalJ (120041) . Nevertheless. 
IRX is a qualitatively different measure of attenuationfrom 
ones that rely on a color excess (e.g. the Balmer decrement 
lKennicut3ll998al and references therein). 

3. ATTENUATION FROM THE IRX-/? RELATION 

The relation between dust attenuation and UV color, the 
IRX-/3 relation (where j3 gives the power-law exponent of 
the UV spectrum, fx ~ A 73 ) is a common tool for the analy- 
sis of the attenuation in galaxies at high-redshift (and low). 
The restframe UV color is more easily measured at high- 
redshift, where it is redshifted into the optical, than many 
other diagnostics (e.g. the Balmer decrement). However, 
at low redshift the calibration has bee n restricted, until re - 
cently, to galaxies with very high SFR (Meurer et al.lll999l) . 
A number of more recent studies have suggested that the re- 
lation becomes significantly scattered when extended to less 
rapidly star-forming galaxies, or even to in dividual OB as- 



sociations, and is shifted to redder UV color ( Bell et al. 2002 
Seibert et al.ll2005l:ICortese et alj2006tlGil de Paz et alj2007 
Boissier et al.l2007l:lDale et al.l2006l) . lKong et al.ld2004l) sug- 
gest, on the basis of stellar population synthesis modeling, 
that the increased scatter and the shift to r edder UV color, 
may be due to the effect of SFH on (5. iBelll (2002) argue that 
these effects are due to radiativ e trans fer through non-trivial 
dust geometries. iPanuzzo et al.1 (120071) find that differing at- 
tenuation of stellar populations with different ages can lead to 
significant changes in the location of galaxies in the IRX-/3 
plane. 

3.1. As a Function o/D„(4000) 

We construct an analog of the IRX-/3 relation for our sam- 
ple galaxies. We use 0A (f -ri) color as a pro xy for /3, since 
we lack UV spectra - see Kong et al. (2004) for a compari- 
son of (3 derived from GALEX colors and from spectra. In 
the left panel of Figure Q] we show the relation between IRX 
and 01 (/-«) color for different ranges of D„(4000). We find 
that the over all scatter in th e IRX- /3 relation is much reduced 
compared to iJohnson et al.l (1200 6). This is due to the much 
smaller errors in UV color with the > 10 times deeper UV 
data presented here. However, there is still no clear trend of 
the scatter in the IRX-/? relation with D„(4000) (our chosen 
SFH indicator) for D„(4000) < 1.7. Such an effect is con- 
strained by these data to be quite small for low D„(4000), 
although it may also be masked by a strong dependence 
of the con version between IRX and the true attenuation on 
D„(4000) (IJohnson et al.ll2007bl). 

As seen in IJohnson et al.1 d2006), for D„(4000)> 1.7 the 
scatter in the relation is greatly increased, and the galax- 
ies have in gener al much redder UV co lor than they do for 
D„(4000) < 1.7. iGil de Paz etall (120071) see a similar scat- 
ter in very UV red galaxies. There is no strong dependence 
of the behaviour at D„(4000)> 1.7 or D„(4000)< 1.7 on the 
SFH indicator being used (e.g. the Ha equivalent width, see 
IJohnson et al.l d2007bl) for examples of additional SFH indica- 
tors available for this sample). 

We also show in Figure [T] the expectation for the relation 
between A/rx a nd 01 (f-n) color for a variety of dust mod- 
els draw n from IWitt & Gordon! ((2000). See Johns onet al.1 
(2007a) for a detailed description of these models and their 
implementation in these color-color diagrams. Note that no 
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FIG. 1. — The IRX-/3 relation shown in two projections. Left: IRX vs. 01 (f — n) color. The color of each point indicates its D„(4000) quintile from lowest 
(purple) to highest (red). Error bars at lower left show the median errors for the entire sample (black) as wel l as the lowest and hi ghest D„(4000) quintiles 
(purple and red, respectively). The black symbols and dashed lines show the effect of dust attenuation from Witt & Gordon 12000), with a SHELL global 
geometry and for combinations of a MW (opensymbols) or SMC (filledsymbols) extinction law and a clumpy (squares) or homogeneous (circles) local dust 
distribution. Symbols run from tv = 0.1 to 10 in steps of 0. 5. Right: Dn(4000) vs. > ( f — n) color. The color of each point indicates its IRX quintile form 
lowest (purple) to highest (red). Overplotted in black are four Bruzual & Chariot (2003) population synthesis models with Z/Zq = 0.4 (open), 2.5 (filled) and 
exponentially declining SFR with time constants tjf =0.1 (circle), 1 (diamond). The symbols are placed at intervals of 500 Myr, from 0.5 to 13 Gyr. The dotted 
lines show Eqs|2]and[3] 



single model matches the data clearly - in particular models 
including a Milky Way like extinction curve are strongly dis- 
favored. The 2 175 A bump in such extinction laws, thought 
to be due to abso rption by ploycyc lic aromatic hydrocarbon 
(PAH) molecules dDraine & Lill2007l) . leads to very little red- 
dening of the UV color with IRX, and perhaps even bluer 
colors with increasing attenuation. Note that small varia- 
tions between galaxies of the particular star/dust geometry, as 
well as the dust extinction law, can easily lead to significant 
scatter in the IRX-/3 diagram. One must also consider here 
the possible effects of a dust attenuation optical depth that 
varie s within the galaxy a s a function of the stellar popula- 
tion dCharlot & Fall 2000), which might include minimizing 
the effect of the 2 175 A bump on the co lor excess as well a s 
significant offsets in the IRX-/3 plane (Panuzzo et alJl2007h . 
Finally, despite the deeper UV data there are still consider- 
able measurement errors in the UV color, and also measure- 
ment a nd systematic errors affecting the calculated IR lumi- 
nosity dJohnson et al.ll2007ah . that can easily lead to consider- 
able scatter. 

The right panel of Figure Q] shows that for a smoothly 
declining exponential star formati on model (SFR~ f~^ TSF ) 
with time constant t$f = 1 Gyr ( Br uzual & CharloJ 120031 : 
iJohnson et alj200 7a). the UV color does not appear to change 
significantly even as D„(4000) increases, except for low- 
metallicity models at late times and high D„(4000). A cor- 
relation between D„(4000) and unattenuated UV color at 
early times would seem to require the presence of significant, 
rapidly declining bursts (e. g. the tsf = 0-1 Gv 1 ' model in Fig. 
[[])■ Indeed, the models of IKong etail d2004l) that show the 
largest deviations from the starburst relation between IRX and 
are those containing significant, short bursts of star forma- 
tion - such SFHs may not be well-represented in this sample. 



The right panel of Figure Q] also shows that the UV col- 
ors of many of the galaxies with red 01 (/-n) color and 
large D„(4000) are consistent with being produced entirely 
by evolved stars from old, passively evolving systems (see 
also lDonase t al. 2007]). As we show injjohnson eildl d2007bl) 
there may also be large contributions to the IR luminosity 
from dust heated by old, redder stars. Thus, the meaning of 
IRX in these galaxies is unclear. For comparison with other 
work we define two cuts, meant to remove such galaxies for 
which IRX is likely be a poor measure of attenuation because 
of the contribution of old stars to both the UV and IR flux. 
These are 

Dn(4000)< 1.95-0.125 x 01 (/-n) (2) 
D n (4000)< 1.625-0.833 x 01 (/-n) (3) 

3.2. Comparison with Other Studies 

A comparison with th e IRX-/3 relat i ons o f IMeurer et al. 
( 1999). lKong et al.l d2004.ISeibert et all d2005l) . ICortese et al. 
(2006), and lGil de Paz et a l. (2007) is possible using the con- 
v ersions from f- n co lor to and from IR X to A f uv given 
in lKong et all d2004l) and lSeibert et alJ d2005l) . We collect the 
various relations here for clarity: 

A /BV =4.37 + 1.74^ (4) 
A fuv = 3.79 + 1.74/3* (5) 
A /BV =5.25 +2.125& (6) 

where 0m e is measured from IUE spectra, s b = 2.286(/- 
n) - 2.096 dSeibert et alJ l2005h . k = 2.201(/-n)- 1.804 
(Kong et al. 2004), we assume 0iue = 0k, and we use the con- 
version from A f m , to I RX given in equati on Q] as was done 
by these authors. The IKong et all d2004l) and IMeurer etail 
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FIG. 2.— IRX vs. ao (/-n) color (an analog of the I RX-/3 relation). Th e 
X'-corrections to z = are determined fol lowing jBlanton & Roweisj (2007ft. 
The IRX-/3 relations of lKong et all <200l) (dottedline) ISeibert et alj <2005h 
(shifted redw a rd by 0.2 mag; solidline), Meurer et al. i 199 91) (dashedline) 
Cortese et al. (2006) (Eq[7J dot - dashedline). and IGil de Paz et"al] <2007l) 
(Eq[8] triple — dot — dashedline) are overplotted (see text for details). Only 
those galaxies that satisfy eqffjare shown. Galaxies that do not satisfy eq|5] 
are shown as open circles. The colors of the points are as in Fig. \U with 
D„(4000) quintiles determined from the entire sample (see Figure 0. Er- 
ror bars in the lower right show the median error for the lowest quintile of 
D„(4000) (purple) and second highest quintile (orange), as well as the me- 
dian error for all galaxies shown (black). 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 

D n (4000) 



FIG. 3. — Top: Deviation in IRX from that predicted by the relation of 
Kongetal. (2004), given the measured (f—n) color, shown as a function 
of D„(4000). Filled circles are galaxies that satisfy eq. [5] open circles are 
galaxies that satisfy eq. [2] but not eq. [3] and crosses are galaxies that do 
not satisfy eq. [2] Bottom: As in the top panel, but expressed in terms of 
magnitudes of attenuation. The dashed lines show the predicted offsets as a 
function of D„(4000) from Kong et al. (2004) (using the the approximations 
given in eqs. 8 and 9 of that work) for 0,0 (f — n) = 0.2 (bottomcurve) and 
°'°{f—n) =1.5 (topcurve). Error bars at top left show the median errors for 
the sample. 

(1999) relations are based on a samp le of starbursting g alax- 
ies observed with IUE and IRAS. The lSeibert et all d2005l) re- 
lation is derived from a sample of galaxies of a wider range 
of types, observed with GALEX and IR AS. An important d if- 
ference between this study and that of Seibe rt et alj d2005l) is 
that the UV color calibration of GALEX has c hanged by ^0.2 
mag su ch that galaxies now have larger f—n, (Mo rrissev et al.l 
(2007), M. Seibert, 2007, private comm unication). We thus 
show the relation of Seibert et al. (2005) shifted by 0.2 mag 



to the red; there is a small effect (< 0.1) of these different cal- 
ibrations on IRX, w hich we ignore. In ad diti on to the three re- 
lations given above. lCortese et aT] d2006l) and lGil de Paz et all 
(120071) give relations that are linear between IRX and (3 or 
(/— «), derived from a sample of nearby galaxies with a wide 
range of SF properties. These are: 

IRX=1.30 + 0.7/3 t (7) 
IRX = -0.18 + 2.05(/-n) (8) 

These five relations are shown in Figure [2] There is rea- 
sonabl e agreement with the (shifted) relation of Seibe rt et alj 
d2005l) . especially at °°(f-ri) = 0.6 and IRX~ 0.7 where 
many of the galaxies are located. However, it is clear that 
the galaxies in our sample, while still obeying a tight relation 
between IRX and 00 (/-n), have a somewhat steeper slope 
than has b een proposed in the p ast. In particular, in agree- 
ment with lGil de Paz et ail d2007l) . the galaxies with very low 
IRX, usually dwarf galaxies of low luminosity, tend to fall 
well below the typical relations. However, several complica- 
tions, in addition to the calibration issues, must be mentioned 
when comparing this sample with others: 

1. Bolometric corrections. - We omit t he IR bolo- 
metric corr ection (BQr = 1.75) used by ISeibert et alj 
d2005l) and iMeurer et al] d 1999ft since we have cal- 
culated to tal dust luminos ities using the IR template 
SEDs of iDale et all d2001l) inst ead of the definitio n of 
iHelouet alJ d!988ft . Similarly. Ilong et all d2004ft. us- 
ing a formula derived from the IDale & Hel ou (2002) 
template SEDs, report IR luminosities ~ 50% higher 
than o btained using the prescription of iHelou etaT] 
( 198 8b, and t hey fi nd a steeper IRX-/? relation than 
IMeurer et alJ d 1999ft for the same sample of gala xies. 
Moreover, as mentioned in iJohnson et al.l d2007al) sig- 
nificant variations in the L(24/zm) to Lj us t ratios are 
possible in the models of IDale et alJ ([2001), and these 
add to the systematic uncertainty in IRX, especially 
given that the mid- to far-IR colors of the sample galax- 
ies do not closely match the temp lates used for the IR 
bolometric correction (Fig. 1 of lJohnson et al.l l2007a). 
Finally, the IRX that we use includes any contribu- 
tion to IR from dust heated by old stars - in particu- 
lar light absorbed in t he optical as opposed to the UV 
dJohnson et al.ll2007bT) - which can cause IRX to over- 
estimate the UV attenuation. 

2. ^-corrections.- A furt her complicat i on is that the 
IMeurer et all d 1999ft and ISeibert etaT] d2005ft relations 
are defined at z = and we must account for the ef- 
fect of the ^-correction to z = 0.1 on the f — n col- 
ors, and to a lesser extent on IRX (due to the depen- 
dence of IRX on /. The extrapolations required to 
go from z = 0.1 (the median redshift of our sample) 
to z = 0.0 can be nearly as large as the size of the 
scatter in the relation: ( 00 (/-n)- 01 (/-«)) = 0.13 
mag with a la dispersion of 0.17 mag, and individual 
/^-corrections from the observed frame f—n color to 

(f—ri) c an be larger than 0. 5 ma g. Note also that the 
method of Bla nton & Roweisl d2007ft relies on spectral 
synth esis modeling of the very type that we are test- 
ing. iBurgarella et al.1 d2005l) have suggested, on the ba- 
sis of GALEX grism spectra and photometry, that the 
observed frame f — n color is a poor tracer of f3 for 
galaxies with z > 0. 1 - we find no significant difference 
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FIG. 4. — SEDs normalized by stellar mass and averaged in bins of several galaxy parameters. From left to right these parameters are IRX, D„(4000), and M*. 
The bins used for color coding are shown above each panel. Error bars give the dispersion at each wavelength band of logM*/i;L„ within each parameter bin - 
the distribution of M t /uL„ is nearly log-normal for the parameter bins chosen here. For vL v and M, in solar units the constant on the )>-axis is 0.99, as for the 
values in Table|2] 



in the distribution of galaxies in Figure [2] if we further 
restrict the sample to z < 0.1 

Not e also that the iMeurer et all (1 19991) and iKong et al.1 
(2004) relations may have been affected by aperture effects in 
the IUE data (M Seibert 2007, private communication), and 
that the GALEX color zero-point error is ~ 0.1 mag. These 
circumstances highlight the difficulty in using the IRX-/3 re- 
lation at high redshift, since it is poorly defined even at low 
redshift and subject to numerous systematic effects, as well as 
random measurement errors that are very important given the 
steepness of the relation. 

In Figure [3] we show the offset (in the IRX direction) of 
the sample galaxies from the Kon g et all (120041) relation (eq. 
|6]l as a function of D„(4000). We also show the offsets of 
Airx from the relation of [Kong et al.l d2004l) for Af uv . There 
is a trend of thes e offsets with D„(4000), as suggested by 
IKong et all (120041) as a signature of the effect of SFH on 
the UV color, although the trend is weak or non-existent for 
D„(4000)< 1.6. The dashed lines in Figure [3] show the pre- 
dictions of IKong et al.l ((2004) (Eq 8 and 9) for two extreme 
values of (/-n). While there is some agreement with these 
predictions at 1.4 <D„(4000)< 1.8 (modulo a constant off- 
set), the predictions do not simultaneously reproduce the off- 
sets at very low D„(4000). The interpretati on is father com- 
plicated by the poor correspondence of the lKong et al.l (2004) 
relation to the galaxies in this sample, evidenced by a large 
offset in Figure |2] 

4. AVERAGE SPECTRAL ENERGY DISTRIBUTIONS 

We have so far considered only the variation in the U V color 
with IRX. iJohnson et al.l d2007al) consider additional colors 
and their variation with IRX. However, it is possible to con- 
struct UV through IR spectral energy distributions (SEDs) for 
the sample galaxies that simultaneously show all colors, in- 
cluding the UV and IR behavior used to infer the SFR and at- 
tenuation. We construct average SEDs for the sample galaxies 
by first nor malizing the galaxy flux at all wavelengths by M* 
as given by lKauffmann et al. (2003a). We then split the galax- 
ies into bins of different galaxy parameters: IRX, D„(4000), 
and M*. These SEDs are shown in Figure [4] and are given 
in Table [2] Note that we have required a 70fim observation 



(i.e. only Lockman Hole galaxies are included) but we do 
not impose a S/N cut at this wavelength; the 70/im points in 
the SEDs, especially for red-sequence galaxies, should thus 
be treated with significant caution. These SEDs will be useful 
for comparison to high redshift galaxies and for theoretical 
studies. 

All three parameters - IRX, D„(4000), and M* - appear 
correlated with significant change of the SED shapes. The 
scatter in SED shape is lowest when the SEDs are binned by 
D„(4000). The change of SED shape as a function of M* is 
driven by the change in SED shape with IRX and D„(4000)- 
M* is correlated with both of these parameters. IRX and 
D„(4000) are only slightly corr elated with each othe r in the 
galaxies in this sample (Fig. 7 of Johnson et al. 2007a), so the 
change in SED shapes with IRX and D„(4000) indicate that 
both attenuation and SFH serve to influence the SED shape. 
T his is shown for individ ual colors (instead of the whole SED) 
in IJohnson et al.l (l2007al) . Note that there is a dependence of 
the 3.6fim stellar mass to light ratio on D„(4000), but little 
dependence of that ratio on IRX. 

To show the separate effects of IRX, D„(4000), and M* on 
the SED shapes we can bin the large sample of galaxies by all 
three parameters simultaneously (after normalizing by M»). 
The results of this process are shown in Figure[5]where instead 
of using quintiles to define the parameter bins we set defi- 
nite bins of width 0.5 in logM* and for D„(4000)we use bin 
edges D„(4000)=[ 1.0, 1.2, 1.4, 1.6, 1.7, 1.9,2.2]. These bins are 
ar ranged to approximate the color-magnitude diagram (Fig. 2 
of IJohnson et al.ll2007al) . Within each of these stellar mass 
and D„(4000) bins we further separate the galaxy SEDs into 

bins of Airx with bin edges A /sx =[0.0,l,2.0,2.5,3.5,4.5,5.5] 
(shown in Fig. [5] as black, purple, blue, green, orange, and 
red, respectively). Many of the SEDs shown are still averages 
of several galaxies, although at the low stellar mass end there 
may be only one galaxy per bin of stellar mass, D„(4000), and 

Airx- It is clear that for a given narrow range of D„(4000) and 
M* there is a significant change in the SED shape due to vary- 
ing Airx ■ At higher stellar mass there is a larger number of 
significantly attenuated galaxies, and of course galaxies with 
higher D„(4000) also tend to have more stellar mass, account- 
ing for much of the change in the SED shape as a function of 
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FIG. 5. — SEDs for different bins of stellar mass (given as log(M*/M0)) and D„(4000), analogous to a color magnitude diagram. Within each of these 
stellar mass and D„(4000) bins the SEDs are normalized by stellar mass and further binned by Airx- See text for details and the definition of the param- 
eter bins. log(M» /Mq ) increases from left to right, while D„(4000) increases from bottom to top. The color of each average SED encodes the Airx bin 
(black, purple, blue, green, orange, and red, respectively) for the six bins defined by A/kx=[0.0,1,2.0,2.5,3-5,4.5,5.5]. 
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FIG. 6. — Derived A^/A;«x (i.e. the dust attenuation law, i|5) for differ- 
ent ranges of D„(4000) (D„(4000)= [1.1,1.2,1.3,1.4,1.5]; bluetored). The 
dashed line shows a A -0 ' 7 attenuation curve. Galaxies are restricted to 
9.5 < log(M* /Mq ) < 11.5. The curves are shifted slightly in wavelength 
for clarity. 

M* (but see fg). 

5. AN EMPIRICALLY DERIVED DUST ATTENUATION LAW 

The change in average SEDs as a function of IRX shows the 
effects of different amounts of attenuation on the stellar mass 
to light ratio, which increase for shorter wavelength. The form 
of this increase shows, in effect, the average attenuation law 
for galaxies. It is possible to derive this attenuation law from 
these data by making several assumptions. It is also neces- 
sary to simultaneously split the galaxies into narrow ranges of 
D„(4000), since the SFH also has a strong effect on the color 
of galaxies. We assume an optical depth: 

T X = C-G(\) (9) 

where we assume that G(A) is universal (the shape of the at- 
tenuation curve) for a given D„(4000) and C provides the nor- 
malization of the attenuation curve. These are separable under 
the additional assumption of reasonably simple geometries. 
In principle this assumption can be tested by repeating the 
analysis detailed below for different ranges of Airx- Then the 
transmitted flux is 

fx=f e- T > (10) 
where f is the intrinsic, unattenuated flux. So, 

A A =2.51og(/ // A ) 
= 1.086 C G( A) 

and the attenuation at a given wavelength may then be written 
in terms of the attenuation at another wavelength as A a, = 
A Al (G(A 2 )/G(Ai)). 

The color of a galaxy, cx u \ 2 = m\ 1 — m\ 2 where m\ is the 
magnitude at wavelength A, is given by the intrinsic color c 
plus the difference in the effective attenuation at each wavel- 
ngth. 

c\ u \ 2 =c + (A Xl -A\ 2 ) 

= c +A Al [l-G(A 2 )/G(Ai)] 

We make the (large) assumption that Airx=A\ 1 = Ai3yo- As- 
suming that c is relatively constant between galaxies for nar- 
row ranges of D„(4000), the slope of the relation between the 
0A (f-m\ 2 ) color and Airx, for a narrow range of D„(4000), 



then give s (1 - G(A2)/G ( 1390)) . This relation is effectively 
shown in Johnson et al. (2007a|) for several different colors 
and ranges of D„(4000), although IRX was used instead of 
A/ R x and many of the colors used 01 n instead of A f. We split 
the sample into bins of D„(4000)(=[1.0,1.2,l. 3,1.4,1.5]) and 
determine the sl opes for each colo r using the ordinary least 
squares bisector (llsobe et al.ll 19901) . In Figure [6] we show the 
resulting derived G(A2)/G(1390) for the different ranges of 
D„(4000). Errors are determined from the range in the slopes 
derived by fitting Airx given 01 (f—tn\ 2 ) and vice-versa. The 
resulting attenuation curves are consistent with a A -0,7 attenu- 
ation law. 

There are a number of caveats to this method forjhe deter- 
mination of the dust attenuation law. Primarily, Airx is not 
the true A 1390, especially for galaxies with redder intrinsic 
spectra. This is likely the reason for the differing shapes of 
the attenuation curves for different ranges of D„(4000). This 
analysis does not exploit the unique advantage of IRX as a 
dust indicator, i.e. that the IR emission constrains the total ab- 
sorbed flux. Indeed, the assumption that A/rx= A1390 is incon- 
sistent with the fact that the derived attenuation laws are non- 
zero at longer wavelengths, which means that galaxies with 
different D n (4000) a re likely to have different r\ (Eq Q] and 
iJohnson et al.l d2007bl) ). Second, G(A) may not be separable 
from C (i.e. the shape of the attenuation curve may be depen- 
dent on the total amount of atte nuation; this can happ en due 
to geometric effects as shown by lWitt & Gordonl (120001) '). An- 
other possible systematic error is that the unreddeded spectra 
may be determined by more than just D„(4000) (e.g. metal- 
licity, contamination due to aperture effects) - c„ may not be 
constant even in small ranges of D„(4000), especially for col- 
ors based on the flux far from the 4000A break. 

Another caveat to the method described above is that G(A) 
may not be universal, for example because of a variation in 
the extinction law between galaxies. Such a variation might 
be expected if the composition of the dust causing the atten- 
uation varies between galaxies. Dust composition may vary 
with metallicity, which is correlated with stellar mass. In par- 
ticular the 2 175 A bump in the Milky Way extinction curve, 
which is not seen in the LMC and SM C extinction curves , 
is thought to be due to PAH molecules dDraine & Lill2007l) 
- these molecules have been shown to be underabundant in 
low-metallicity dwarf galaxies, although the reason for this 
is still unclear. A variation in the attenuation law with mass 
would be of special importance for interpreting the location 
of gal axies in the specific SFR - stellar mass plane in terms of 
SFHs (lLabbe et al.l2007l;lNoeske et al.ll2007l) and f or compar- 
ison of UV to Ha derived SFRs (ISalim et al. ,2006). For these 
reasons we consider the attenuation laws derived as above but 
for bins of D„(4000) and stellar mass (i.e. cells in Figure 0. 
These are shown in Figure|7] The smaller number of galaxies 
in each bin results in a weaker constraint on G(A2)/G(1390)) 
than in Figure [6] No significant trends with stellar mass are 
seen, although the uncertainties are large. A larger sample of 
galaxies is required to accurately test this hypothesis and, in 
addition, split the sample into, e.g., bins of Airx to explore the 
coupling of attenuation law shape to the amount of attenuation 
in galaxies. 

6. ATTENUATION AS A FUNCTION OF METALLICITY AND 
SURFACE BRIGHTNESS 

At the local level the amount of dust attenuation is due to 
the column density of dust, r = J adustndmt^ where adust is 
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FIG. 7.— As in Figure|6] the derived A x /A IR x for different ranges of D„(4000) (D„(4000)= [1.1, 1.2, 1.3, 1.4, 1.5]); bluetored) and stellar mass. Each 
panel shows a different stellar mass range. We only show those stellar mass and D„(4000) ranges for which there are > 10 galaxies. The dashed line shows a 
A -0 7 attenuation curve. The curves are shifted slightly in wavelength for clarity. 

star formation. Note also that AGN identified in the optical 
spectra on the basis of emission line ratios are excluded from 
this s ample with metallicity measurements (Tremonti et al. 
2004). The left and center panels show that there is a strong 
correlation between IRX and E* even when the metallic- 
ity is held nearly constant; metallicity and surface mass den- 
sity ar e also correlated wit h each other, as shown in the right 
panel. ICortese et alj d2006l) found a similar trend of IRX with 
H-band surface brightness (a close proxy for stellar mass sur- 
face density), although they did not simultaneously control for 
metallicity. The correlation of IRX with metallicity is very 
weak, for constant surface mass density. However, there is an 
overall correlation of IRX with metallicity induced primarily 
by the correlation of metallicity with surface mass density. 
The lines in Figure [8] show those expected from equation [Til 
for the quintiles of the parameter used for the color coding 
(i.e., points of a certain color should fall between two adjacent 
lines). We have assumed f g = 0. 1 and require 77 = 0.7 to fit the 
data reasonably - this relation is very sensitive to the choice of 
-q. The lines in the middle panel, for different quintiles of the 
metallicity distribution, show that the effect of the metallicity 
variations is expected to be weaker than the surface mass den- 
sity variations, since the range of metallicity is only 1 order of 
magnitude while the range in surface mass density is 2 orders 
of magnitude. Note that the shallower-than-expected slope of 
the data in this panel may well be due to the anti-correlation 



the absorption cross section, nd us t is the dust volume den- 
sity, and t is the unit pat h length along the line of sight (e.g. 
Wang & Heckman 1996). Simplistically, one might write 
ndusd = ^gasfdust where fdua is the global dust to gas ratio 
and Egns is the column density of gas along the line of sight 
to the star or cluster of stars. We assume that fd„ s r oc Z, 
where Z is the g as-phase metalli c ity. W e obtain Z from the 
logO/H values of iTremonti et alj (|2004) via the solar values 
(logO/H) Q = 8.69 and Z Q = 0.02, assuming that the oxygen 
abundance is a tracer of the total metal abundance within a 
factor of 2 or 3. In semi-analytic modelling it is typical to 
make the assumption fd us t oc Z, and furthe rmore to use the 
globa l values of gas densi ty and metallicity ( Som erville et alJ 
2001; Mar tin et all [2006). This approach involves several 
other assumptions: that the effect of stars being embedded 
within the gas can be neglected, that spatial variations in both 
gas surface density and metallicity are such that their product 
averages to the global value, and that all stars (of all ages) 
are affected equally. Such a relation between surface density, 
metallicity, and attenuation has been investigat ed empirically, 
using glob al quantities of resolv ed galaxies, by lBoissier et alJ 
d2007l) and lCortese et al.1 d2006l) . 

We consider here the relation of the globally measured at- 
tenuation to the global metallicity and surface density for the 
sample galaxies. We do not have measurements of the gas 
surface density, and use instead the stell ar mass surface den- 
sity. Following Kauf fmann et al.1 (l2003bl) th is is calculated as 
E iiZ _= 0.5M*/(7r/?5 ) where M* is given by Kau ffmann et alJ 
(2003a) and R50 is the radius (in kpc) en closing 50% of the 
Petr osian light in the z b and. iBelll d2003l) . Appendix B, (see 
also ICalzetti et al.ll2007l Appendix A.l) makes a quantitative 
estimate of the relation in this case, obtaining: 



T/„ v ,1550 = l-77?Z/^/(l-/g)S» 



(11) 



where r/ MV ,i550 is the optical depth at 1550A, Z is the metal- 
licity, f g is the gas fraction, E* is th e ste l lar m ass surface 
density, and r\ is a factor of order unity dBelll ([2003) determine 
77 = 0.7 to fit the attenuation-luminosity correlation) to account 
for, e.g., relative star dust geometry. Using the approximation 
Tf U v = t / „,, 1 55o(1390/1550)" ' 7 and equation 1, r/„ v ,i550 can be 
related to IRX. 

In Figure [8] we show three projections of the relation be- 
tween IRX, gas-phase metallicity, and E*. z . The number of 
galaxies for which this is possible is small (219), due to the 
lack of metallicity measurements - these galaxies have strong 
emission lines and typically have large ratios of recent to past 



of f g , which we have assumed to be constant, with surface 
density. 

Indeed, it is possible to estimate E P from the star form ation 
rate density, assuming the relation of lKennicuttl (fl998b) holds 
globally for these galaxies. Gas densities are calculated via 

= lO 6 [(lO 4 /5)SFR g /(7r ^ o )]-'- 4 M -kpg- 2 (12) 

using the median SFR derived by Brinchm ann et al.1 d2004l) . 
and /?5o as above (the conclusions are unchanged if we in- 
stead use the u band half light radius as an estimate of the 
extent of star formation). The mean of the ratio Eg^/E*^ 
is approximately 0.16 (hence the assumption f g = 0.1 made 
above). There is a significant trend of this ratio with surface 
mass density E*^, such that galaxies with larger E* )Z have 
lower E ?J t/E*. z - E^ spans a much smaller range than E*. z . 
This explains the steeper slope predicted by eauation fTTl 

The se results may be compared to those of Boissier et al. 
d2007|), who found a cor r elation of IRX with metallicity (see 
also lBoissier et aUl2004t iHeckman et~ai1ll998l) . These stud- 
ies did not consider the simultaneous effect of the surface 
density-metallicity correlation, and a trend with metallicity 
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FIG. 8. — Three projections of t he relation between IR X, gas-phase metallicity, and surface mass density for the 219 galaxies of the sample for which metallicity 
measurements are available from Tremonti et al. 12004). In each large panel the color coding is given b y the quintiles of the distribution shown above that panel. 
In the left panel the black lines show (bottomtotop) the relation of IRX to metallicity expected from eq. 1 1 1 I for the values of S*.,- that define, in increasing order, 
the limits of the bins used for the color coding. In the middle panel the lines are similar but here the metallicity is fixed at the limits of the bins used for the color 
coding. 



may thus be due to this additional component of the optical 
depth. What is perhaps more puzzling is that Boissi er ;t alj 
d2007l) did not find a good correlation between th e local IRX 
and th e lo cal HI(+H2) su rface density (but see Bua t & Xul 
(1996) and IXu et alj dl997h ). The relation between gas sur- 
face density and stellar surface density is highly scattered and 
poorly constrained, although we have assumed that the global 
values are simply related by f g . If the stellar surface mass 
density is better correlated with IRX than the gas surface den- 
sity we must ask why. Is the stellar mass density tracing the 
history of metal and dust production, such that higher densi- 
ties naturally lead to higher attenuations? This would seem 
to be accounted for by considering narrow ranges of metal- 
licity. One component of the analysis that we have treated 
only superficially is that stars are necessary to measure IRX 
and hence the attenuation - a region of large gas and dust col- 
umn will not contribute to the measured attenuation if there 
are few stars, and conversely a large number of stars in a low 
dust column region will cause a low measured IRX, even if 
the global gas and dust surface density is large. IRX in fact 
measures the sum of the local Ti gas Z weighted by the local 
stellar luminosity. 

This relation between surface-density, metallicity, and IRX 
is very closely related to the relations between IRX and lumi- 
nosity and between IRX and stellar mass - these, in turn, are 
key to interpreting the color-magnitude diagram o f galaxies in 
terms of models for the b uild-up of stellar mass (Lab be et al.l 
20071 lNoeske"eTal]l2007h . If the sizes of late-type galaxies 
are nearly constant or correlated with the total SFR, then the 
gas surface densities are also correlated with total SFR. In this 
case, the relation between IRX and gas surface density, a more 
direct probe of the physical conditions related to attenuation, 
is almost trivially re lated to the well-k nown correlation be- 
tween IRX and SFR (Ma rtin et aT1 l2005). A similar argument 
may be made for the stellar mass surface density and total 
stellar mass. Interestingly, the scatter in the relation presented 
here does not appear to be significantly less th an the scatter in 
the IRX-L.s F relation (see lJohnson et al.ll2007bl) . 

If, however, the measured attenuation does prove a use- 
ful probe of the global S ?flJ Z, then the measured distribu- 
tion of attenuation may provide a strong constraint on mod- 
els of galaxy evolution, since both gas surface density and 
metallicity are intimately related to star birth and star death 



dMartin et al.ll2 006). 



7. CONCLUSIONS 



1. For D„(4000) < 1.7 we find little evidence that scatter 
in the IRX-/3 diagram is caused primarily by variations 
in UV color with the SFH of the stellar population, as 
suggested by Kong et al. (2004). However, IRX may be 
affected by the SFH, which could mask such an effect. 

2. We assess the utility of the IRX-/3 relation for determin- 
ing dust attenuation at high redshift using UV colors 
alone. While the scatter is small, the large slope makes 
the derived IRX very sensitive to even small errors in 
the UV color (~ 0.1 mag). These may be random (due 
to measurement errors) or systematic, for example due 
to the method used to /^-correct the data. Furthermore, 
significant uncertainty in the determination of the bolo- 
metric IR luminosity makes even the low redshift rela- 
tion uncertain. 

3. We use the UV through IR data to show how the en- 
tire spectral energy distribution (not just a single color) 
varies as a function of the relevant galaxy parame- 
ters IRX, D„(4000), and M*. This results in high- 
quality, low resolution 1375 A through 70/im average 
SEDs normalized by for different ranges of IRX, 
D„(4000), and M*. Such average SEDs show a vari- 
ation of the 3.6/im mass to light ratio with D„(4000), 
when s tellar masses are e s timated from optical data fol- 
lowing |K^^mannetan d2003a|). They also show the 
strong variation in the UV and IR specific luminosities 
with these quantities, although the variation with M* is 
driven largely by the correlation of SFH and attenua- 
tion with M* (which might be expected to change with 
redshift). 

4. We have used these average SEDs to derive a low res- 
olution dust attenuation curve for blue-sequence galax- 
ies, for several different ranges of D„(4000). The ac- 
curacy of the attenuation curves is limited by the vari- 
ation of the relation between IRX and the true attenua- 
tion A f m , but the derived attenuation curves are consis- 
tent with a A -0,7 law. To investigate possible changes in 
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the attenuation law as a function of metallicity or PAH 
abundance we repeat the analysis, splitting galaxies by 
their stellar mass. We find no significant differences as 
a function of mass although the uncertainties are large. 

5. IRX is correlated with both stellar mass surface 
density and gas-phase m etallicit y (as determ i ned by 
iKauffmann et all (l2003bl) and iTremonti et all (120041) . 
respectively), although the latter correlation appears 
weak for a given narrow range of stellar mass sur- 
face density. To the extent that the stellar mass sur- 
face density is related to the gas surface density, IRX 
then probes the causes of dust attenuation and this re- 
lation will allow the accurate calculation of attenuation 
in models of galaxy formation. This, in turn, will al- 
low for the consistent comparison of photometric data 
to these models. Conversely, the attenuation can be 
used to probe the variation in the surface density and 
metallicity as a function of time through observations 
of high-redshift galaxies in the restframe UV and Mid 
to Far-IR. 
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TABLE 2 

Average" Broadband SEDs Normalized by M, in Bins of IRX, D„(4000), and M* 



Bin < {og(uF„/M«)+d° > (CT[logOF„/M*)]) 
Min. Max. | ol n 0A u ol g 0A r 01 ; A z 01 3.6^m 4.5^m 5.8^m 7.8^m 24/^m 70/im 

Binned by IRX 


-0.46 0.61 
0.61 0.87 
0.87 1.09 
1.09 1.39 
1.39 2.63 


0.19 (0.467) 0.18 (0.395) 0.30(0.201) 0.66(0.115) 0.77(0.071) 0.81 (0.061) 0.78 (0.042) 0.06(0.035) -0.22(0.035) -0.32(0.074) -0.03 (0.197) -0.59 (0.369) 0.76(0.415) 
-0.44 (0.456) -0.39 (0.382) -0.09 (0.264) 0.39 (0.062) 0.56 (0.034) 0.62 (0.029) 0.64 (0.021) -0.03 (0.027) -0.28 (0.030) -0.46 (0.058) -0.21 (0.272) -0.80 (0.404) 0.46 (0.468) 
-0.71 (0.551) -0.57(0.424) -0.11 (0.113) 0.35(0.050) 0.54(0.029) 0.60(0.024) 0.62(0.019) -0.02(0.026) -0.27(0.033) -0.44(0.069) -0.23(0.349) -0.83(0.514) 0.49(0.549) 
-0.95(0.414) -0.76(0.307) -0.23(0.087) 0.27(0.043) 0.48(0.026) 0.55(0.023) 0.58(0.017) -0.04(0.021) -0.28(0.027) -0.48(0.053) -0.30(0.289) -0.79(0.376) 0.52(0.529) 
-1.17(0.300) -0.89(0.194) -0.32(0.061) 0.19(0.033) 0.41 (0.023) 0.49(0.021) 0.52(0.017) -0.03(0.024) -0.27(0.036) -0.42(0.058) -0.13(0.242) -0.54(0.292) 0.77(0.366) 


Binned by D„(4000) 


1.06 1.31 
1.31 1.45 
1.45 1.66 
1.66 1.85 
1.85 2.20 


0.37 (0.201) 0.39 (0.158) 0.46(0.095) 0.77 (0.068) 0.84(0.048) 0.88 (0.042) 0.85 (0.030) 0.20(0.020) -0.06(0.022) -0.09(0.032) 0.43 (0.056) 0.01 (0.080) 1.38(0.083) 
-0.18(0.206) -0.12(0.139) 0.08(0.066) 0.45(0.040) 0.59(0.027) 0.65(0.022) 0.65 (0.017) 0.03(0.013) -0.21 (0.017) -0.31 (0.021) 0.18(0.041) -0.29(0.080) 0.97(0.116) 
-0.63(0.262) -0.50(0.151) -0.16(0.063) 0.30(0.034) 0.49(0.025) 0.56(0.020) 0.58(0.016) -0.05(0.010) -0.29(0.016) -0.46(0.020) -0.12(0.054) -0.62(0.094) 0.69(0.132) 
-1.09(0.291) -0.89(0.158) -0.35(0.071) 0.20(0.020) 0.43(0.016) 0.51 (0.013) 0.55 (0.011) -0.10(0.010) -0.35(0.015) -0.59(0.015) -0.52(0.085) -1.10(0.140) 0.17(0.294) 
-1.52(0.118) -1.28(0.078) -0.47(0.103) 0.13(0.015) 0.39(0.012) 0.47(0.011) 0.52(0.009) -0.14(0.014) -0.40(0.020) -0.67(0.019) -0.86(0.054) -1.51 (0.120) -0.17(0.255) 


Binned logM* 


0.00 9.70 
9.70 10.50 
10.50 10.75 
10.75 11.00 
11.00 12.50 


0.54(0.557) 0.55 (0.356) 0.59(0.317) 0.93 (0.087) 0.99(0.054) 1.02 (0.045) 0.97 (0.032) 0.20(0.033) -0.11 (0.041) -0.10(0.070) 0.13 (0.165) -0.29 (0.277) 1.13 (0.280) 
-0.29(0.667) -0.16(0.416) 0.14(0.126) 0.53 (0.069) 0.67 (0.039) 0.72 (0.032) 0.71 (0.024) 0.06(0.026) -0.22(0.033) -0.30(0.060) 0.06(0.257) -0.43 (0.325) 0.86(0.488) 
-0.58(0.446) -0.46(0.309) -0.10(0.086) 0.34(0.040) 0.53(0.021) 0.59(0.017) 0.61 (0.013) -0.01 (0.016) -0.27(0.024) -0.43(0.045) -0.11 (0.237) -0.66(0.329) 0.65(0.406) 
-0.84 (0.385) -0.69 (0.284) -0.22 (0.068) 0.27 (0.030) 0.47 (0.019) 0.54 (0.016) 0.57 (0.013) -0.05 (0.019) -0.28 (0.029) -0.49 (0.040) -0.28 (0.244) -0.79 (0.365) 0.50 (0.415) 
-1.12(0.300) -0.97(0.220) -0.40(0.122) 0.15 (0.023) 0.39(0.015) 0.46(0.013) 0.50(0.010) -0.13(0.013) -0.33(0.026) -0.59(0.033) -0.49(0.194) -1.08(0.347) 0.23(0.366) 



Because vF u /M* is close lo log-normally distributed for the parameter bins given here, the mean of the log of this quantity is an appropriate statistic to report, as opposed to the mean or the log of the mean. C = 0.99 



